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On the Optimality of Napoleon Triangles 
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Abstract. An elementary geometric construction known as Napoleon’s theorem produces an 
equilateral triangle built on the sides of any initial triangle: the centroids of each equilateral 
triangle meeting the original sides, all outward or all inward, comprise the vertices of the 
new equilateral triangle. In this note we observe that two Napoleon iterations yield triangles 
with useful optimality properties. Two inner transformations result in a (degenerate) triangle 
whose vertices coincide at the original centroid. Two outer transformations yield an equilateral 
triangle whose vertices are closest to the original in the sense of minimizing the sum of the 
three squared distances. 

1. INTRODUCTION. In elementary geometry, one way of constructing an equilat¬ 
eral triangle from any given triangle is as follows: in a plane the centroids of equilat¬ 
eral triangles erected, either all externally or all internally, on the sides of the given 
triangle form an equilateral triangle, illustrated in Figure 1 [I]. This result is generally 
referred to as Napoleon’s theorem, notwithstanding its dubious origins — see [2] for 
a detailed history of the theorem. We will refer to these constructions as the outer and 
inner Napoleon transformations and the associated equilateral triangles as the outer 
and inner Napoleon triangles of the original triangle, respectively. Conversely, given 
its outer and inner Napoleon triangles in position (i.e. they are oppositely oriented and 
have the same centroid), the original triangle is uniquely determined [3]. A fascinating 
application of Napoleon triangles is the planar tessellation used hy Escher: a plane can 
he tiled using congruent copies of the hexagon defined hy the vertices of any triangle 
and its outer Napoleon triangle, known as Escher’s theorem [4]. 






Figure 1. An illustration of (left) the Fermat point F, outer Torricelli configuration and outer 

Napoleon triangle Aj^n^n^n, and (right) inner Tonicelli configuration Aaj,BtCt inner Napoleon 
triangle of a triangle Aj^gc- Note that centroids of Torricelli configurations, Napoleon 

triangles and the original triangle all coincide, i.e. c{Aabc) = ^(A^TgT^T) = c{Aj^n^nqn) = 
cCA^ySyCr) — ^(^Aj^gffCff)- 
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Equilaterals built on the sides of a triangle make a variety of appearances in the clas¬ 
sical literature. Torricelli uses this construction to locate Fermat’s point minimizing the 
sum of distances to the vertices of a given triangle, now called the Fermat-Torricelli 
problem [5]. The unique solution of this problem is known as the Fermat-Torricelli 
point of the given triangle, located as follows [6]. If an internal angle of the triangle is 
greater than 120°, then the Fermat point is at that obtuse vertex. Otherwise, the three 
lines joining opposite vertices of the original triangle and externally erected triangles 
are concurrent, and they intersect at the Fermat point, see Figure 1 . The triangle defined 
by the new vertices of the erected equilateral triangles is referred to as the Torricelli 
configuration [7, 8]. 

In this paper we demonstrate some remarkable, but not immediately obvious, op¬ 
timality properties of twice iterated Napoleon triangles. First, two composed inner 
Napoleon transformations of a triangle collapse the original to a point located at its 
centroid which, by definition, minimizes the sum of squared distances to the vertices of 
the given triangle. Surprisingly, two composed outer Napoleon transformations yield 
an equilateral triangle optimally aligned with the original by virtue of minimizing the 
sum of squared distances between the paired vertices. More precisely, for any triangle 
^ABC with the vertices ' A, S, C G we will say that the triangle Aa'b'C is an 
optimally aligned equilateral triangle of Aabc if it solves the following constrained 
optimization problem: 


minimize ||A — A'\f -j- \\B — B'\f -f \\C — C'\^ 
subject to II A' - = P' - C'f = \\B' - C 


where A' ,B' ,C' G and ||.|| denotes the standard Euclidean norm on As we 

show below, this optimization problem has a unique solution so long as A, B, C are 
not collinear. 


2. TORRICELLI AND NAPOLEON TRANSFORMATIONS. For any ordered 
triple X = [xi,X2,X3] G R^'^, let Rx denote the rotation matrix corresponding to a 
counter-clockwise rotation by 7r/2 in the plane, defined by orthonormal vectors n and 
t, in which the triangle Ax formed by x is positively oriented (i.e. its vertices in 
counter-clockwise order follow the sequence ...—?-l—)-2—^3—^1 —t'...), ^ 


Rx : = [ n, t ] 


0 -1 

1 0 


[ n, t f, 


( 2 ) 


where ^ 

_^^2, _ JG {z G S'^-^ln^z = 0}, if X is collinear 

“ 1otherwise, ' “ |p(n) , otherwise, 

where P(n) : = 1^ — nn^ is the projection onto (the tangent space of at 

point n G and I^^ is the d x d identity matrix. Note that Ax is both positively 

and negatively oriented if x is collinear. Consequently, to define a plane confaining 

'Here, R denotes the set of real numbers, and M'* is the d-dimensional Euclidean space, 
denotes the transpose of matrix A. 

^For any trivial triangle Ax all of whose vertices are located at the same point we fix Rx = 0 by setting 

TiAi- = 0 whenever x = 0. 

I|x|| 








such X we select an arbitrary vector t perpendicular to n in (3). It is also convenient 
to have c(x) denote the centroid of Ax, i.e. c(x) ; = ^ X]i=i 

In general, the Torricelli and Napoleon transformations of three points in Euclidean 
d-space can be defined based on their original planar definitions in a 2-dimensional 
subspace of containing x. That is to say, for any x G select a 2-dimensional 
subspace of containing x, and then construct the erected triangles on the side of Ax 
in this subspace to obtain the Torricelli and Napoleon transformations of x. Accord¬ 
ingly, let Tj_ : —)• and Nj_ : —)• R^'^ denote the Torricelli and Napoleon 

transformations where the sign, + and —, determines the type of the transformation, 
inner and outer, respectively. One can write closed-form expressions of the Torricelli 
and Napoleon transformations as: 

Lemma 1. The Torricelli and Napoleon transformations of any triple x G R^*^ on a 
plane containing x are, respectively, given by 


where 


T±(x) = (^iK± 

N±(x) = i(Kx + T±(x)) , 




X, 


'o 

1 

l' 



0 

-1 

l' 

1 

0 

1 

® Id 

and L = 

1 

0 

-1 

1 

1 

0 



-1 

1 

0 


(4) 

(5) 


( 6 ) 


Proof One can locate the new vertex of an equilateral triangle, inwardly or outwardly, 
constructed on one side of Ax in the plane containing x using different geometric 
properties of equilateral triangles. We find it convenient to use the perpendicular bi¬ 
sector of the corresponding side of Ax, the line passing through its midpoint and being 
perpendicular to it, such that the new vertex is on this bisector and a proper distance 
away from the side of Ax. 

For instance, let y = [yi, y2, Ys]^ = T^(x). Consider the side of Ax joining Xi 
and X 2 , using the bisector, bi 2 : = 4(xi + X 2 ), to locate the new vertex, ya, of in¬ 
wardly erected triangle on this side as 


Ys = bi2 + —Rx(x2 - Xi) , (7) 

where Rx (2) defines a counter-clockwise rotation by ^ in the plane where x is posi¬ 
tively oriented. Note that the height of an equilateral triangle from any side is times 
its side length. Hence, by symmetry, one can conclude (4). 

Given a Torricelli configuration y = [yi, y 2 , Ya] = (x), by definition, the ver¬ 

tices of associated Napoleon triangle z = [zi, Z 2 , Za]^ = Nj_(x) are given by 

11 1 

Zi = ^(Yi +X 2 +Xa), Z 2 = -(xi +y 2 + Xa) and Za = -(xi +X 2 +ya), (8) 

which is equal to (5), and so the result follows. ■ 


^Here, (g) denotes the Kronecker product [9]. 







Note that the Torricelli and Napoleon transformations of x are unique if and only if 
X G is non-collinear. If, contrarily, x is collinear, then is both positively and 
negatively oriented and for d > 3 there is more than one 2-dimensional suhspace of 
containing x. 

Remark 1 ([3]). For any x = [xi, X 2 , Xa]^ € R^*^, the centroid of the Torricelli con¬ 
figuration y = [yi,y 2 ,y 3 ] = Tj_(x), the Napoleon configuration z = Nj_(x) and 

the original triple x all coincide, 

c(x) = c(y) = c(z), (9) 

and distances between the associated elements of x and y are all the same, i.e. for any 
f G{1,2,3} 


l|y*-x ,||2 = l|yj-Xj-lla- (10) 

An observation key to all further results is that Napoleon transformations of equi¬ 
lateral triangles are very simple. 

Lemma 2. The inner Napoleon transformation of any triple x = [xi, X 2 , X 3 ]"'" G 
R^'^ comprising the vertices of an equilateral triangle Ax collapses it to the trivial 
triangle all of whose vertices are located at its centroid c(x), 

N+(x) = I 3 ( 8 )c(x), (11) 

whereas the outer Napoleon transformation N_ reflects the vertices 0 /Ax with respect 
to its centroid c(x), ^ 


N_(x) = 2 • I 3 (g) c(x) — X. (12) 

Proof Observe that the inwardly erected triangle on any side of an equilateral triangle 
is equal to the equilateral triangle itself, i.e. (x) = x, and so, by definition, one has 
(11). Alternatively, using (5), one can obtain 

N+(x) = ^(Kx-f T^(x)) = i(Kx-f x) = l 3 ( 8 )c(x), (13) 

where K is defined as in ( 6 ). 

Now consider oufwardly erecfed equilaferal triangles on the sides of an equilateral 
triangle, and let y = [yi,y 2 )y 3 ]^ = T_(x). Note that each erected triangle has a 
common side with the original triangle. Since Ax is equilateral, observe that the mid¬ 
point of the unshared vertices of an erected triangle and the original triangle is equal to 
the midpoint of their common sides, i.e. |(yi -|- Xi) = |(x 2 + X 3 ) and so on. Hence, 
we have T_ (x) = Kx — x. Thus, one can verify the result using (5) as 

N_ (x) = i (Kx -f T_ (x)) = -^ (Kx -f Kx - x) = 2 • I 3 (g) c(x) - x. (14) 

00 H 


^Here, I 3 is the R® column vector of all ones, and ■ denotes the standard array product. 



Since the Napoleon transformation of any triangle results in an equilateral triangle, 
motivated from Lemma 2, we now consider the iterations of the Napoleon transfor¬ 
mation. For any A: > 0 let —)• denote the k-th Napoleon transformation 

defined to be 

N^+^=N±oN^, (15) 

where we set N^ : = id, and id : —)• is the identity map on 

It is evident from Lemma 2 that: 

Lemma 3. For any x G and k>l, 

N^+^(x) = I 3 (g) c(x), and N^+^(x) = N^(x). (16) 

As a result, the basis of iterations of the Napoleon transformations consists of Nj_ and 
N^, whose explicit forms, except N^, are given above. Using (5) and (12), the closed- 
form expression of the double outer Napolean transformation can be obtained as: 

Lemma 4. An arbitrary triple, x = [xi,X 2 ,X 3 ]"^ G R^'^ gives rise to the double outer 
Napoleon triangle, : R^'^ —R^'^, according to the formula 

Nt(x) = jX+1t^(x). (17) 

Proof. By Napoleon’s theorem, N_(x) is an equilateral triangle. Using (5) and 
Lemma 2, one can obtain the result as follows: 

N^(x) = N_(N_(x)) =2-l3(g)c(x)-N_(x) = 2-l3(g)c(x)-i(Kx-hT_(x)),(18) 

o 

= ^ (Kx-f x) - ^ (Kx-f T_ (x)) = ^ (Kx-T_ (x)) = ^x-f iT+(x) ,(19) 

where K is defined as in (6). ■ 

Notice that N^ (x) is a convex combination of x and (x), see Figure 2. 


3. OPTIMALITY OF NAPOLEON TRANSFORMATIONS. To best of our 
knowledge, the Napoleon transformation N^ is mostly recognized as being a function 
into the space of equilateral triangles. In addition to this inherited property, has 
an optimality property that is not immediately obvious. Although the double inner 
Napoleon transformation is not really that interesting to work with, it gives a 
hint about the optimality of : for any given triangle yields a trivial triangle all 
of whose vertices are located at the centroid of the given triangle which, by defini¬ 
tion, minimizes the sum of squared distances to the vertices of the original triangle. 
Surprisingly, one has a similar optimality property for : 

Theorem I. The double outer Napoleon transformation N^(x) (17) yields the equi¬ 
lateral triangle most closely aligned with Ax in the sense that it minimizes the to¬ 
tal sum of scjuared distances between corresponding vertices. That is to say, for any 




Figure 2. (left) Outer, Aj^n^nqn , and double outer, , Napoleon transformations of a triangle 

^ABC- (right) The double outer Napoleon triangle Aj^ogDQD is a convex combination of the original 
triangle Aabc and its inner Torricelli configuration ■ 


X = 


[xi,X2,X3]^ 


G N^(x) is an optimal solution of the following problem 


minimize 


3 





subject to ||yi - yaf = ||yi - ya 


I|y2 - ya 


( 20 ) 


where y = [yi, y 2 j ys]"*" S R^'^- Further, ifn. is non-collinear, then (20) has a unique 
solution. 

Proof Using the method of Lagrange multipliers [10], we first show that an optimal 
solution of (20) lies in the plane containing the triangle Ax- Then, to show the result, 
we solve (20) using a proper parametrization of equilateral triangles in R^. 

The Lagrangian formulation of (20) minimizes 


3 

^(y,Ai,A2) = ||x* -y*||2 +Ai(||yi -yalla - ||yi -yaliz) 

i=l 

+ Aa^llyi-yalla - ||y2-yaliz) , (21) 

where Ai, Aa G R are Lagrange multipliers. Necessary conditions for the locally opti¬ 
mal solutions of (20) is ® 


y^{y 1 Aij A 2 ) — 2 


(yi - Xi) + Ai(y3 - ya) + A2(yi - y2) 

(y2 - xa) + Ai(y2 - yi) + A 2 (y 3 - yi) 

(ya - X3) - Ai(y3 - yi) - A2(y3 - y2) 


= 0 , ( 22 ) 


from which one can conclude that an optimal solution of (20) lies in the plane con¬ 
taining Ax- Accordingly, without loss of generality, suppose that Ax is a positively 
oriented triangle in R^, i.e. its vertices are in counter-clockwise order in R^. 

In general, an equilateral triangle Ay in R^ with vertices y = [yi,y2, ya]"*" € R® 
can he uniquely parametrized using two of its vertices, say yi and ya, and a binary 


'’Here, Vy denotes the gradient taken with respect to the coordinate y. 







variable k G { —1,+1} specifying the orientation of Ay; for instance. A: = +1 if 
Ay is positively oriented, and so on. Consequently, the remaining vertex, y 3 , can be 
located as 

1 \/3 

¥3 = ^(ji+y2) + k—'R^/2iy2-yi), (23) 

where 'R.t ^/2 = 

Hence, one can rewrite the optimization problem (20) in term of new parameters as 
an unconstrained optimization problem as: for yi, 72 G and k G {—1,1}, 


is the rotation matrix defining a rotation by 7r/2. 


minimize 


1^1 - yiL + 11^2 - 72112 + ||x3 - Myi - M^yal 


(24) 


where M : = |I — A:-^R,r/ 2 , and I is the 2x2 identity matrix. Note that M + 
= I, M^M = MM^ = I and = -M^. 

For a fixed k G { — 1,1}, (24) is a convex opfimizafion problem of yi and 72 be¬ 
cause every norm on M" is convex and compositions of convex functions wifh affine 
fransformafions preserve convexify [11]. Hence, a global optimal solution of (24) oc¬ 
curs where fhe gradienf of fhe objective function is zero af. 


(I -f M^M) 

(M 2 )^ 

7 i 


xi -f M^xa' 

M 2 

(I -h MMT)_ 

/2. 


_ X2 + Mxa _ 


which simplifies fo 


21 

-M 

7 i’ 


Xi + M^Xg' 

-M^ 

21 

72. 


X2 + Mxg 


(25) 


(26) 


Nofe fhaf fhe objective function, /(y), is sfrongly convex since ifs Hessian, V^/(y), 
satisfies 


vV(y) = 


21 -M 

-M^ 21 




(27) 


which means thaf for a fixed k G { — 1, +1} the optimal solution of (24) is unique. 
Now observe that 


1 
3 

hence the solution of linear equation (26) is 


21 

M 

21 

-M' 


I 0 

MT 

21 

-MT 

21 


0 I 


yi 

1 

21 M 

Xi+M^xa 

1 

2xi+2M'^X3 +Mx 2+M^Xg 

.^2. 

“3 

MT 21 

_ X 2 +MX 3 ^ 

“3 

.MTxi+(M2)% + 2X2 + 2Mx3_ 


2xi + + Mx2 

2X2 + M'^Xi + Mx3 


j(^2+ 3 

Here, substituting y^ and 72 back into (23) yields 

Xi -f X2 + X3 


2 (Xi + ’'l +X3 ) ^ R^/ 2 (x3 _ X2 ) 

^il±^)+fe^R./2(x3-X3) 


73 = ^ ( X3 + 


+ -Xi) 


(28) 

(29) 

(30) 


3 


(31) 








































Thus, overall, we have 


y = 



V3, 


-Kx+A: 

2 2 


R-^/2) Lx 


f|x + iT_^(x), if fc = +1, 
||x+iT_(x),iffc = -1. 


where K and L are defined as in (6). Recall that Ax is assumed to he positively ori¬ 
ented, i.e. Rx = R 7 r/ 2 , and so it is convenient to have the results in terms of Torricelli 
transformations Tj_ (4). As a result, the difference of y and x is simply given hy 


i(T^_(x) - x), if A; = +1, 
|(T_(x) - x) , if A; = -1 


Finally, one can easily verify that the optimum value of k is equal to +1 since the 
distance of x to its inner Torricelli configuration T_(_ (x) is always less than or equal to 
its distance to the outer Torricelli configuration T_(x). Here, the equality only holds 
if X is collinear. Thus, an optimal solution of (20) coincides with the double outer 
Napoleon transformation, N?_(x) (17), and it is the unique solution of (20) if x is 
non-collinear. ■ 


As a final remark we would like to note that our particular interest in the optimal¬ 
ity of Napoleon triangles comes from our research on coordinated robot navigation, 
where a group of robots require to interchange their (structural) adjacencies through 
a minimum cost configuration determined by the double outer Napoleon transforma¬ 
tion [ 12 ]. 
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